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Abstract. We compute the characteristic Cartan connection associated with a system of 



third order ODEs. Our connection is different from Tanaka normal one, but still is uniquely 
associated with the system of third order ODEs. This allows us to find all fundamental 
£f~) , invariants of a system of third order ODEs and, in particular, determine when a system 

of third order ODEs is trivializable. As application differential invariants of equations on 
circles in K" are computed. 
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1 Introduction 

> 

*0 1 1.1 Differential equation as a structure on a filtered manifold 

\Q ■ 

The main purpose of this article is to study geometry of systems of ordinary differential equa- 
tions of third order. The geometry of ordinary differential equations or, more generally, of 
\ differential equation of finite type is based on the general theory of geometric structures on 

filtered manifolds. First it was developed by Tanaka in [§J [TO]. Recall that a filtered manifold 
is a smooth manifold M equipped with a filtration of the tangent bundle TM compatible with 
the Lie bracket of vector fields. At any point x € M the associated graded vector space gr T X M 
can be endowed with a Lie algebra structure. This nilpotent Lie algebra m is called a symbol of 
a filtered manifold (at the point x). In the paper we consider only the so-called filtered manifolds 
of constant type, assuming that the graded nilpotent Lie algebras gr T X M are isomorphic to each 
other for all points x 6 M. 

By a symbol of a geometric structure on M we understand a graded Lie algebra g with 
the negative part g_ = Xw<o 0* which is equal to the symbol m of the filtered manifold M 
of constant type. Here the Lie algebra m is the subalgebra of a so-called universal Tanaka 
prolongation g(m). Roughly speaking, this means that g(m) is the maximum among graded Lie 
algebras which satisfy the condition "for any element X € Qi, i > the equality [X, g_] = 
implies X = 0" . 

An arbitrary equation £ can be viewed as a surface in jet space. The canonical restriction of 
the contact distribution on jet space defines the structure of filtered manifold on £. 



1.2 The problem of equivalence 



One of the main problems in the theory of differential equations is the problem of equivalence. 
Two differential equations are called equivalent if one can be transformed to another by a certain 
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change of variables. We consider equations up to point transformations, i.e. we allow arbitrary 
changes of both dependent and independent variables. 

First classical approach to the equivalence problem of ODEs was developed by Sophus Lie. 
In [5] he obtains partial results about second order ODEs. The complete answer was given later 
by Tresse [11]. Invariants of the third order ODEs were computed by Chern in his paper pQ. 
A modern approach to the equivalence problem of ODEs can be found in the papers [2] and [3], 
where characteristic Cartan connection was constructed for the one equation of arbitrary order 
and for the system of ODEs of the second order. 

The general approach to the equivalence problem for the holonomic differential equations 
can be find in [3]. The key fact there is the existence of a full functor from the category 
of holonomic differential equations to the category of Cartan connections. This reduces the 
equivalence problem for differential equations to the equivalence problem for the corresponding 
Cartan connections. 

1.3 Normalization of Cartan connections 

Let P be the principal H -bundle. Let w be a Cartan connection of type (G,H), where G is 
a Lie group with a semisimple graded Lie algebra q and H is a parabolic subgroup of G with 
the Lie algebra f). In the paper [10] Tanaka built a set of normal Cartan connections on the 
principal bundle P as follows. He used the scalar product defined with the help of the Killing 
form to construct adjoint Lie algebra codifferential d*. Then a Cartan connection is normal iff 
the structure function C : P —> Hom(A 2 g_,0) belongs to the kernel of the operator d* and the 
structure function has not negative components. As usual define a Laplacian A = d*d + dd* . 
The structure function C decomposes as C = H{C) + A(C). The component H(C) is called the 
harmonic part of the structure function. The key fact about it is that H(C) is the fundamental 
system of invariants (see Definition [5] for details). In the case of the geometry of holonomic 
differential equations the Lie algebra g is not necessarily semisimple. However in [3] is shown 
that we still can find the scalar product on g such that the normal Tanaka conditions define the 
unique Cartan connection associated to a holonomic differential equation. 

In this paper we associate with every system of ODEs of third order a characteristic Cartan 
connection which differ from a normal Tanaka Cartan connection. The reason for doing this is 
a relation between conformal geometry and geometry of the system of the third order ODEs. 
Conformal manifold is determined by the family of conformal circles, which was shown by 
Yano [12]. Each conformal circle is determined by the point on it, the direction and the curvature, 
i.e. by the point in the third jet space. The system of appropriate differential equations of 
the third order gives us the bridge between the conformal geometry and the geometry of the 
differential equation. It is appeared that a characteristic Cartan connection, which is built in 
the paper, is in close relations with the normal conformal Cartan connection. The relation of 
the conformal geometry and the geometry of third order ODEs is the topic of the next paper. 

The paper is organized as follows. In Section [2] we naturally associate the system of the 
third order ODEs with the pair of distributions. This pair of distributions give rise to the 
filtered manifold associated with the system of the third order ODEs. We write down the 
symbol of the system of ODEs of the third order, the notion of adopted coframe and adopted 
Cartan connection. The problem of equivalence is considered in Section [3l When we working 
in the case of semisimple Lie algebras and normal Cartan connections, the harmonic part of the 
curvature gives us the fundamental system of differential invariants. We show that in general 
case fundamental differential invariants are contained in the Ker<9 part of the curvature, where 
d is the Lie algebra cohomology differential. In Section [J] we build the characteristic Cartan 
connection uniquely associated to the the system of ODEs of the third order. This connection 
allows us to obtain the results about equivalence third order equations and to describe the 
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structure of the fundamental invariants of the system of third order ODEs. In particular, this 
answers the question "When is the given system trivializable?" explicitly. 



2 Geometry of the systems of third order ODEs 

Consider an arbitrary system of m ordinary differential equations of third order: 

y'i'( x ) = fi{y'j(x),y' k ( x ),yi(x),x), (1) 

where i, j, k, I = 1, . . . , m and m > 2. 

We associate a filtered manifold with this system in the following way. Let J 3 (M m+1 ,l) 
be the third jet space of unparametrized curves. Then the equations ([TJ can be considered as 
a submanifold £ in J 3 (M' m+1 , 1). We introduce the following coordinate system on the surface £ : 



{x,yi,-. -,y m ,pi =y'i,--- ,p m = y'm><ii = y'L ■••iQm = y'm) 



There is a natural one-dimensional distribution E whose integral curves are the lifts of solu- 
tions of equations (pQ). Let 7rf be the canonical projection from the surface £ to the first jet space 
J 1 (M m+1 , 1). We denote a kernel of a differential d^ as V. In coordinates distributions E, V 
have the form: 

zp / ' 9 9 9 ri 9 \ r, I 9 

\dx dyi dpi dqi I \c% 
where i, j = 1, . . . , m. 

Define a distribution C as the direct sum of the distributions E and V . Then C and its 
subsequent brackets define a filtration of a tangent bundle T£: 

C = C~ l C C~ 2 C C~ 3 = T£, 

where C"^ 1 = C' { + [C~\ C' 1 ]. 

It is easy to see that the symbol of the filtrated manifold £ is a nilpotent Lie algebra m 
isomorphic to the Lie algebra of vector fields generated by 

Id d d d 

1 \ dx ^ dyj ^ dpj ' dqi 

Let Auto(m) be a subgroup of grading preserving elements of the group Aut(m). The elements 
of the group Auto(m) which preserve the splitting E © V form subgroup Gq. So the splitting 
E © V of the distribution C defines Go-structure of type m. The action of the group Go on m 
is completely determined by its action on m_i. The latter has the following form in the basis 



dx 



+ Pi dyj + 1jd Pj i dq % } : 

a£R* Be GL r 



a 
B 

The symbol g is the universal Tanaka prolongation of the pair (m, go)- It has the following 
form: 

= (s/ 2 (R) x 5 l m (R))A(V 2 ®W). 

In other words, g is equal to the semidirect product of the Lie algebra 0/2 0&) x 0^ m (^) an d an 
Abelian ideal V. The ideal V has the form V2 <S> W, where V2 is an irreducible s/2-module of 
dimension 3 and W = R m is the standard representation of g[ m (R). 
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Let us fix a basis of the Lie algebra 5I2 and sl2-module V2. Let x, y, h be the standard basis 
of an algebra 5I2 with relations: 

[x, y] = h, [h, x] = 2x, [h, y] = -2y. 

This basis can be represented in the following way: 

x= (o J)' *=(J -1)' o)- 

Let uo, ui, V2 be a basis of the module V2 such that x.vi = v\, x.v\ = vo, x.vq = 0. 
Define the grading of the Lie algebra Q as follows: 

01 = (y), 00 = (h,Qt m ), 0_i = (x) + (« 2 ® W), 

Q_2 = {V1®W), Q~3=(V ®W). 

To build a natural Cartan geometry associated to the equation ([I]) we will use the fact [7j 
that under some additional conditions (which are satisfied for geometric structures arising from 
holonomic differential equations, see [3]) there exists a full functor from the category of Go- 
structures of type m to the category of Cartan connections of type (G,H), where G and H are 
the Lie groups with Lie algebras q and f) respectively which are determined from Go in natural 
manner. The group G is a semisimple product: 

G = (SL 2 (R) x GL m (R)) A (V 2 <g> W) , 

and the group H is the following subgroup of G: 

H=(^ jOxA, oeK*, beR, A e GL m (R). 

Note that the corresponding subalgebra f) is exactly the nonnegative part of the Lie algebra q: 

i>0 

Definition 1. We say that a coframe {w!_ 3 ,o;^ 2 ,a;l 1 ,a; a .} on 8 is adapted to equation ([1]) if: 

• the annihilator of forms V: 

• the annihilator of forms u) l _ 3 , u l _ 2 , ^-i is E\ 

• the annihilator of forms is G -2 . 

Let W: P — > £ be a principle ff-bundle and let cJ be and arbitrary Cartan connection of type 
(G, H) on P. Connection U can be written as: 

oj = + uf_ 2 v i ® e « + w l _it>2 ® ej + cJ^x + ZJ^/i + uJj-ef + 

Definition 2. We say that a Cartan connection u on a principal if-bundle 7f is adapted to 
equations ([I]) , if for any local section s of W the set 

S Ld x ,S U)_ 1} 8 W_ 2 ,S W_ 3 ) 

is an adapted co- frame on £. 

We have described the set of Cartan connection adapted to the system of third order ODEs. 
However, we can chose the representative in different ways. The next two sections are devoted 
to the building of a canonical connection which we call characteristic. 
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3 Characteristic Cartan connection 

and fundamental differential invariants 

As in Section [2] let W: P — >• £ be a principle .ff-bundle and let U be and arbitrary Cartan 
connection of type (G, H) on P: 

oj = cJ*_ 3 fo <8> + u) l _2 v l ® e i + ^-i v 2 ® ej + w x x + ZJ^/i + + uJ y y. 

Let f2 = d uJ + 2 [w, uJ] be the curvature of the Cartan connection ZJ: 

f2 = fi?_ 3 vo <g> ej + H!_ 2 ui ® e i + ^-1^2 <8> ej + Q x x + O^/i + + Vt y y. 
Definition 3. The structure function of a Cartan connection w is a function 

C : P -> Horn (A 2 g_,g), 
which is defined by 

C{p){g\,92) = %{^p l {gi)^p 1 {g2))- 

We can obtain the structure function of a Cartan connection explicitly. Let {ei, . . . , e n +fc} be 
a basis of Lie algebra g such that {e n +i, . . . , e n+ k\ form a basis of the subalgebra h. In our case 
{e n+ i, . . . , e n+ fc} = {/i, y, e^}. An arbitrary element <^ 6 Hom(A 2 g_, g) defined by constants C^, 
where 

tp(ei,ej) = ^C^-efe, l<i,j<n. 

k=l 

The structure function C : P — >■ Hom(A 2 g_,g) defines functions C^(p). If 

then the functions C^(p) can be found from the decomposition of the curvature tensor Q, in 
terms of forms u^: 

Let be one of the 2-forms fil_3, ^_2> ^-i> ^x, ^h, fi!-- We can write it explicitly as: 

3 3 
p,q=l p=l 

Then Cl 1 ^ ,JjJ^_ ] and Q l \O x , oJ- p ] are the coefficients of the structure function of the Cartan 
connection uj. The grading of Lie algebra g induces degree of the coefficients Q l [lJ :> _ g ,uj k _ p \ and 

Definition 4. We say that Cartan connection associated with the equation ([I]) is characteristic 
if the following conditions on a curvature is satisfied: 

• all coefficients of degree < 1 are equal to 0; 

• in degree 2 we have fi/j[57 x A57!_x] = 0, [u^Ao;^] = 0, £l x [^xA5j!_ 2 ] = 0, r2„ 1 [cJ x .AtJ*„2] = 0; 
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• in degree 3 we have £l y \uJ x A oJ'l 1 ] = 0, Qhfex A ^-2] = 0' \px A w* 2 ] = 0; 

• in degree 4 we have £l y \uj x A cJL 2 ] — 0. 

In other worlds these conditions define the subspace U and Cartan connection is characteristic 
if and only if it belongs to U. 

Theorem 1. There exists a unique characteristic Cartan connection associated with the equa- 
tion (H|). 

Proof. We will proceed with parametric computations of characteristic Cartan connection in 
the forth section of the paper. We will fix a section s : £ — > P and prove that locally for 
every equation there exists a unique Cartan connection u with structure function pullback 
s*C : £ Hom(A 2 g_,g) takes values in the space U. Now we show that the characteristic 
Cartan connection is uniquely globally defined with this data. 

Take a covering U a of the space E and construct a Cartan connection u a on each trivial 
fibre bundle Tr a : U a x H — > U a . Let s a and sp be the trivial sections of the fibre bundles 7r a 
and 7rg. Let ui a = s* a oj a and ujp = spuip. Since forms oj a and ojp are uniquely defined there 
exists a unique function 

(fa/3 :U a nUf3 ^ H, 

such that 

Up = Ad (^})w a + <p* a pUH, 

where ujjj is Maurer-Cartan form of the Lie group H. The functions ip a p uniquely define 
a principle i?-bundle with the Cartan connection u. 

In order to prove that the structure function C of the Cartan connection u) takes values in 
the space U it is sufficient to show that U is Ad(i/)-invariant. 

Note that the action of Go preserves the zero condition on the structure function of the 
characteristic connection. We need only to check that the space U is exp(y)-invariant or equally 
ad(y) invariant. The action of the element y has degree one. Conditions on the curvature of 
the Theorem [1] are ad(y)-invariant up degree 2, since all components of degree less than 2 are 
equal to zero. Finally, the conditions of degree 3 and 4 are ad(y)-invariant, since the coefficients 
0, y \uJ x AuP_i\, Qfifex AtJ^K ^j[w x AaJ^L 2 ] anc ^ ^y^x^ 1 --^ can be- obtained only from f^[aJ x .AZJl_i], 
J7 x [lJ x . AuJ l _ 2 ]i Qj\px Aujkj], Oft[cJ x Aa7!_ 2 ] an d ^j/[^ x A uf_i\ which all are zero for characteristic 
Cartan connection. This ends the proof of a global existence of the form co. ■ 

Let V be an arbitrary finite-dimensional vector space and let / be a smooth function 
/ : P — > V. Denote by Lo(f) the space of all functions of the form (f,v*), where v* € V* 
and by L(f) the algebra generated by elements from Lo(f) and all their covariant derivatives. 
For example, the algebra L(C), where C is structure function of the Cartan connection u), 
consists of local invariants of the connection u. 

Definition 5. We say that functions /j are the fundamental system of differential invariants 
for the structure with Cartan connection uj if L(fi) = L(C). 

The key to calculation of the fundamental system of differential invariants is to determine 
which parts of the curvature are expressed through another. In [3] it is shown that fundamental 
invariants of holonomic differential equation lie in non-negative harmonic part of the curvature 
of the normal Cartan connection. In general we have approximately the same situation: there 
is one to one correspondence between fundamental differential invariants of the characteristic 
Cartan connection and H^(q^,q) part of the structure function. Here iJl(g_,g) is the non- 
negative part of the second Lie algebra cohomology group. 
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Proposition 1. Let uj be a Cartan connection of type (G,H) on a principal H -bundle P, where 
(G, H) is an arbitrary pair of Lie group and its subgroup. Assume that the Lie algebra g is 
a graded Lie algebra of the Lie group G with the negative part g_. Assume that a structure 
function of oj takes values in subspace W C Hom(A 2 g_,g) and has only components of posi- 
tive degree. Then a Ker d n W part of the structure function forms a system of fundamental 
differential invariants. 

Proof. The algebra of differential invariants is generated by the structure function coefficients. 
We will use the Bianchi identity to show that some coefficients of the characteristic Cartan 
connection curvature are obtained from the image of the operator d. 

Let e, be the basis of the Lie algebra g, Xi be the corresponding fundamental vector fields 
on P and u l be the dual coframe. We can write the Cartan connection oj in the form: 

uj = u l ei. 



Assume that the Lie algebra g has structure constants A\-. That means that: 
[ej, ej] = A^je^. 



Write the curvature of the Cartan connection oj in coordinates: 

n = Cgw i AoJe k . (2) 
Then the following equality is fulfilled: 

duj k = (Cg-Aj)a/ Aoj j . 
Now apply the Bianchi identity df] = [fi,u>] to the equation (J2j) : 

(dC- \ 
-Q^-Ul A u)i A u)je k + d Ui A ujje k + C^uji A d u)je k J = [e k , e^Ui A ojj A uj t . 

Express the covariant derivative of the structure function: 

Q(JP, 

-q^-Ui A 0Ji A uijCp = (-C^duJk Aui - C v k] uj k A dwi + C^A^Ui A Uj A uii)e p 
= Cf,(C7§ - A!j)ui A uj A uie p + C^A p M uji A Uj A uie p . 

We get that: 

Q(JP 

-^jrOJi A uii A UjCp - ClftijLOi A ujj A w/e p 

= C^AijOJi A LOj A oj t e p + CijA^LOi A Uj A LOie p . (3) 

If we take the Hom(A 3 g_,g) part of (|3|) (i.e. assume that uji G gl) we get that the right side of 
the ([3]) is exactly the Lie cohomology differential. 

On the right side of (J3|) coefficients have the same degree as in the curvature. On the 
other hand coefficients on the left side have an increased degree. So, we have obtained that 
coefficients which are mapped to the im d can be expressed through the covariant derivative of 
the coefficients of the lower degree. This proves the proposition. ■ 

Remark 1. Note that if intersection of W and Im d is zero then subspace Ker dD W is generated 
by representatives of Hi. (g_,g). 
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Theorem 2. The following invariants are fundamental differential invariants for the system of 
third order ODEs: 



wo; 



tr 

df 



dp _ d_df_ 
dpi dx dqi 

1 dp df k 



1 df 1 df 



3 dq k dqi 



ih) 



+ 



dyi 3 dq k dpi 
5 d fdp\ df k 



1 d dp 1 d 2 dp 

2 dx dpi 6 dx 2 dqi 



2 

27 



tr 

df 
dq k 



d 2 P 



dqidq k 

3 



18 dx \dq k J dqi 
d d 



d^ 1 



+ 



d Pj dq 3 dq k 



H x 



dqk dx j 



where 



H: 



l 



d 2 P 

6(m + 1) \dq l dq^ 



and 



H 3 



d 

Qgk 



1 

4m 



H, 



-idf 
dqi 



1 5/* d df k 
mdqk 'dx~dq7 



dp 

dp 1 



d dp 1 dp dp 
dx 9(f 3 dq k dq l 



Proof. We will use Proposition [TJ The fundamental differential invariants is in one to one 
correspondence with the cohomology group if+(g_,g). For the case of the system of ODEs of 
the third order the Lie cohomology group q) was studied in [6]. The main result of that 

work is that the space has the following decomposition: 



Degree 


Space 


-1 


1 
2 
2 
3 
4 
3 


ug® A 2 (W*)®W 
vl®Sl(W*)®W 
v° 4 ® A 2 (W*)®W 
v% ® A 2 W* <g> W/V 2 ® W* 
x* ®Ry®sl(W) 
v%®S 2 {W*)®W 
x* ® Ry 2 ® fll(W) 

if m = 2 



Here vjj! is the lowest vector of corresponding (fc + l)-dimensional s/2-module 

Now we list the result table with the corresponding invariant. We start from degree 2 since 
all part of curvature of degree less than 2 is zero. 



Degree 


Space 


Part of the curvature 


Invariant 


2 
2 
3 
4 
3 


x* ®Ry ®sl(W) 
v%®S 2 (W*)®W 
x* ®Ry 2 ®Q\(W) 
v%<E> S 2 (W*) 
v® if to = 2 


^!,K A wi 3 ] 
A o;i 2 ] 


W 2 
/ 2 

/ 4 
= 



Corollary 1. T/ie system (pQ) is equivalent to the trivial one via point transformations if and 
only if all invariants I2, W%, W3, I4 vanish identically. 

Example 1 (Differential equations on circles in R n ). As application of the previous results we 
compute invariants of the system of third order ODEs on circles in Euclidean space. 
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Lemma 1. Let E be the (m + 1)- dimensional Euclidean space with the orthonormal basis 
{eo, • • • , e n } and the coordinates {ro, r±, . . . , r n }. Then the equation of circles in E parametrized 
by the coordinate ro is: 

m 

r'i = 3r\ J ~ , i = l,...,m. 



l+£r 2 



m 

r 

J 



J'=l 

This equation is invariant under conformal transformations of E. 

Proof. Let the curve R(t) = (ro(t), ■ ■ ■ , r n (t)) be a circle. Assume now that ro(t) = t. We have 
R(t) = a(t)R(t) + b(t)R(t), (4) 

since R(t) is 2-dimensional curve. Next, b(t) = in our parametrization, since 
= r (t) = a{t)r {t) + b(t)r {t) = b(t). 

To determine a(t) note that 
(R{t)-C,R(t)-C) = d 

for some constant d and C € E. Differentiating, we get: 

(R(t),R(t)-C) = 0, 
(R(t),R(t)-C) = -(R(t),R(t)), 
(R\t),R(t) -C) + 3(R(t), R(t)) = 0. 

Now substitute ([!]) into previous formula: 

(a(t)R(t) + b(t)R(t), R(t) -C) = -3{R{t) , R(t)) , 
(a(T)R(t),R(t) -C) = -a(t)(R{t),R(t)) = -3(R(t) , R{t)) . 

We get that 

(R(t),R(t)) 

a(t) = 6 : : . 

(R(t),R(t)) 

Substituting a(t) into (jU) we get our equations. ■ 

Proposition 2. For differential equation on conformal circles invariants Wi, I2, W3 vanish 
identically. Invariant I4 has the following form: 

/ T \i 1 f' TjTi 
( J 4) v - = 



'1 2 \ ™ .0 2 / m \- 



Proof. The proof is straightforward applying of the formulas from Theorem [2j ■ 

Remark 2. There are other equations satisfying W2 = I2 = W3 = 0. For example, it is an 
union of a system on circles in R n ~ k and a system of k trivial equations. It would be interesting 
to characterize geometrically the class of such equations. 
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4 Parametric computation 

of the characteristic Cartan connection 

Consider a system of third-order ordinary differential equations of the form 

(yT = f{^y j ,iy k )',(y 1 )"), 

where i,j = l,...,m with m > 2. It determines a holonomic differential equation £ C 
J 3 (]R m+1 , 1). Let us use the following coordinate system on the equation £ : 

x,yi, . . .,y m ,pi = y'i, ■■•,Pm = y' m Ai = y", ■■■,qm = y f m- 

We choose a coframe 9 on the surface £: 
9 X = dx; 

0-1 = dq l - f{x,y,p,q)dx, i = l,...,m; 
9 % _ 2 = dp 1 — (f dx, i = 1, . . . , m; 
9 l _ 3 = dy % — p l dx, i = 1, . . . , m. 

To connect our computation on the surface £ with the principle bundle P let us use the following 
uniquely defined section s: £ — > P with relations: 

s*uJ h = mod (9^,9^,6^), 
s*ul x = -9 x mod <0L 3 ,0- 2 ,0Li>. 

Define a pullback uj : T£ — > g by the formula uj = s*uj. Let fl be a curvature tensor of UJ, and 
let f2 = s*£l. We see that 

O = n l _ 3 v ®ei + SVl 2 vi <g> e» + ^1^2 ® e» + fi^ + O^/i + Of + 0. y y 
= (du> l _3 + oj x f\ uj % _2 + 2u h A u l _ 3 + A uxL^Vo <g) e; 
+ (do;L 2 + c^e A o;Lx + A o;i 2 + 2w y A uj' 1 _ 3 )vi <g) e, 
+ (do/.! - 2o; ft A wlj -1- ujj A ar^ + 2oj y A u;!_ 2 )f2 ® 

+ (dcu x + A lo x )x + (dcj/j + w x A ojy)^ + (dujj + col A u j) e i + (duly ~ A 
An arbitrary Cartan connection adapted to equation (pQ) has the form: 
coU = 9l 3 , 

U l 2 = a ^_ 2 + A)9 j _z, 

■J i -ij(' J , • irf • r;^ ,. 

w x = — 9 X + Dj9~!_2 -^jO—31 

co h = FJ l 9l x + F" V 2 + F j V ;i . 

W ) = Gf ^ + G}- 1 ^! + G)- 2 9 k L 2 + Gp9l 3 , 

Lo y = H X 9 X + Hj 1 9l l + Hj 2 9l 2 + Hr 3 0L 3 . 

In degree of the curvature we have two nonzero components: 

ni 3 mod (£_2 A 9^ 2 , 9- 3 ) = X A 9l 2 - a)9 x A 9 j _ 2 , 
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ni 2 mod (0_ 2 , #_ 3 ) = e x a - p)6 x a &_ v 

Assume these two equalities is zero and get ct- = 5j and /3j = 5j. 

We have three nonzero components in degree 1. The first component is: 

nt 3 mod (9-2 A 9-3, 9- 3 A 9- 3 ) 

= -9 X A A)eL 2 + DjOL 2 A 9L 2 + (r; r (), A 9 3 i 3 + G^T 1 ^ A 9 j _ 3 + 2^7^ A 9l 3 . 

The second component is: 

Qi 2 mod (9- 2 A9- 2 ,9- 3 ) 

= A)9 X A #i 2 + Dj-flij A 0Li — 9 X A B)9 j _ 2 + Gf9 x A 9 3 _ 2 + G i ^ x 9 k _ 1 A 9 j _ 2 . 

The third component is: 

ni x mod(0_ 2 ,0_ 3 ) 

= A ei x + A flij - 2F7 i 9 j _ 1 A 9l x + Gf9 x A 9 3 _ x + G i - k ~ 1 9 k L x A 9 j _ v 



After applying zero conditions to these parts of the curvature we obtain 

Proceed now to the second degree 

mod (9- 2 A9- 2 ,9- 3 ) 

8 P dA i dA l 

= w-Ox A 9 3 + 2—^9 x A 9 J _ 2 + 2—±9K A 9 3 + C% A 9l 2 - 2Fr 2 9 3 A 9\ 
opj ax dqk J J 

+ G^ 2 9 k : 2 A 9 3 _ x + 2H X 9 X A 9i 2 + 2Hj 1 9 3 _ 1 A 9i 2 + G\ X 9 X A B k 9 3 _ 2 . 

We have: 

Of , n dA) 

Assuming the previous tensor is zero, we obtain: 



1 

i Ak 



«Li \0 X A 9 3 _ 2 ] = 7t- — h 2—^ + C\ + 2/f B + 2A\A 
1 " J apj ax J J 



dp dA\ . ,\ 

CI- = - + 2—^ + 21P + 2A\A k . 
3 ydpj dx k 3 J 

Next curvature component contains all second order invariants: 



St_ 2 mod (0_ 2 A 0_ 3 , 0-3 A 0- 



3 



dA i dA i 

= ~^0 X A 9 3 _ 3 + -j-LO 1 - , A £3 — 9 X A C)9 3 _ 3 + ^0i 3 A 0! , + Gf 0* A A{9l, 
+ 2H X 9 X A 9i 3 + 2HJ 1 9l 1 A 9l 3 + G i ^ 2 9 k _ 2 A 9 j _ 2 + G i ^ 1 9~ 1 A Aj9 l __ 3 . 
In coefficient ^!_ 2 [^i A 0i 3 ] we get invariant I 2 

Qi 2 [9 k , A 9l 3 ] = ^ - ^-4 + 2fl^ = ^ + 2ffr^j + 2F- 2 4. 
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Explicitly, the invariant I 2 is the following: 

o 2 r 



tr 



dqjdq k 



where tro is a traceless part of the tensor. 
In the coefficient 

dA i 

ni 2 [e x a ei 3 ] = -cj-^r + a% a j + 2Hk $j 

we obtain a so-called generalized Wilczynski invariant. As shown in 0, a part of differential 
invariants of systems of ODEs comes from its linearisation. As in [2], we call them generalized 
Wilczynski invariants. In our case we have two Wilczynski invariants of degree 2 and 3. We 
denote them as W 2 and W3 respectively. The second degree generalized Wilczynski invariant is 
the following: 

{dp d dp 1 dp df k 
W 2 = tr hr + - 



dpj dx dqj 3 dqk dqj 
Normalizing the trace of previous tensor to zero we obtain: 

H* = -± (?f + 3^ + 3A\A k 
Am \dpi dx 

It remains to compute only sl 2 x Q\ m part of the curvature in degree 2. 

n x mod (6-2 A 9-2, 9 -3) = Ej9 x A 9 j _ 2 + 2Fj0 x A 9 j _ 2 . 
Assuming that it vanishes identically we get the following condition: 

Ej = ~2F7\ 
We have: 

n h mod {9- 2 , #_ 3 ) = Ff% A 9l x - 9 X A ^flj 1 . 

The condition A 6>!_i] = gives equality F~ 2 = HJ 1 . 

Assuming the trace of the tensor Q, l _ 2 [9 :> _ 1 A 9^_ 3 ] is equal to zero we get: 

F -2 _ H -i_ 1 94 
k k 2(m + 1) dq k ' 

The last part of degree 2 calculation is: 



dA 

n} mod (fl_ 2> #_ 3 > = ^ a ^ + G^ 2 e, a 0^. 



We obtain G l - k 2 = from condition Ttj[uj x A w^J = 0. 

Proceed now to the degree 3. The first part of degree 3 we need to compute is £l z _i- 

BP ■ dB i 
mod (9^2 A 9- 3 ,9-3 A 9-3) = ^-9 X A J _ 3 + —^9 k 2 A 9 3 _ 2 

dyi dp k 

Q(ji Q(ji 

+ ^9 X A 0i, + -=-^0* 1 A 0: L, - 2iA~V „ A 0Li - 2Fr V 2 A BL0* 2 
ax c% J ^ 
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+ G jk 3d -3 + G^ 2 e k : 2 a B)ei 2 + a'r'o,. a c 3 k ei 3 + 2Hj 2 el 2 a ei 2 . 



Wilczynski invariant W3 appears as the f^f^ A 9 J _ 3 ] coefficient: 

df dC) . , 

J- + —1 + A\C k + 2H X A). 



Direct computation shows that: 

ni .„ „.• , df ldf l df k d df 2 d 2 df 2 fdf 



dyi 3 dq k dpi dx dpi 3 dx 2 dqi 27 \ dq> 
Adf d df 2d ( df\ df 



2S 1 H' x 

9dq k dxdqi 9 dx \dq k J dqi j 

Denote invariant f2!_i[^i A 6 3 _ 3 ] + \ ^f^2 as W3. Invariant W3 is equivalent to the fundamental 
invariant fi!_i[#3i A 0^3] ■ It means that after replacing A 0^_ 3 ] with W3 the system would 

remain fundamental. Explicitly the Wilczynski invariant W3 is: 

dp 1 df df k 1 d df 1 d 2 df 



9y-' 3 dq k dpi 2 dx dpi 6 dx 2 dqi 

df\ 3 1 df d df 5 d /df \ <9/ A 



27 \dqi J 18 dq k dx dqi 18 dx \dq k J dqi 

An expression (§^j) 3 here is the third power of the matrix Note that invariant W% has 
known analogue in the case of one differential equation of third order: 

df_ \df_df_ _ \±_df_ + I d2 df 2 fdf\ 3 ldfddf 



dy 3 dq dp 2 dx dp 6 dx 2 dq 27 \ dq J 3 dq dx dqi 

The reader can find this invariant for example in Chern work pQ; also see Sato and Yoshikawa [8]. 
Let us compute the third degree normalization conditions. 



F~ 3 e x a ei 2 + —j—e x a el 2 + -^—e^ a e j _ 2 -9 x a Afei 2 . 



n h mod (L 2 AL 2 ,0_3> 

ax agfc 
Thus: 

d\F~ 2 

Normalizing this coefficient to we obtain: 

FT 3 = Hf - d -Sl. 

J J dx 

Next, 

dA i dG % '~ 2 
SI) mod (0_ 2 A 0_ 2 , 6L3) = 2 A 6> x + A 9 k _ 2 + A 9 k _ 2 



+ G*'% A Gjf^L 2 + G l £'e k _ 2 A Gf : 
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We have: 



ft Ai r!C l '~ 2 

JI± 3 , jk . . n i,x n l-2 n i,-2 n l,x 



Assuming this coefficient is equal to we get: 

G * ~^~^r~ Gi G > k + Gik G >- 

Finally, 

dH x dHj 1 
Q y mod (0_ 2 , 0_ 3 ) = -Q—0-1 A ^ + A 0^ 

dH 1 f)fj 
+ -^-^-i a 0ii + Hj l -Le x A 0* x + ff/ 2 ^ A 

The coefficient A fl^] is the following: 

dqj dx k dqj 3 

Normalizing it to we obtain: 

H -2_dH* dH- 1 H ^dfl 
3 dqj dx k dqk ' 

The last coefficient we need in degree 3 is Q y [6 k i 1 A O^}: 

dHr 1 9H -i 

3 - = 0. 

dq k dqj 

In the degree 4 we need to compute only one coefficient of curvature: 

■ „ 9Hr\ h . dH, 



-2 



Q y mod (6. 2 A 0_ 2 , 0- 3 ) = ^—^2 A ^ + I ^0- 2 A 0i x + —?-6 x A 0i 2 

j 



dH~ ftfj 

The Cartan connection coefficient £l y \6 x A 0^2] nas the following form: 
dH* | dHf | ^Q/* 

Assuming it is equal to we get: 

R _ 3 _ dH- | dHf | ^g/fc 

Finally, invariant I4 is the tensor 1 A 9 -2}- 

OH' 1 dHr 2 x j 
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